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Abstract 
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£SJ , We introduce diquarks as separable correlations in the two-quark Green's 

function to facilitate the description of baryons as relativistic three-quark 
■ bound states. These states then emerge as solutions of Bethe-Salpeter equa- 

tions for quarks and diquarks that interact via quark exchange. When solving 
these equations we consider various dressing functions for the free quark and 
diquark propagators that prohibit the existence of corresponding asymptotic 
Q_i| states and thus effectively parameterize confinement. We study the impli- 

O 1 cations of qualitatively different dressing functions on the model predictions 

for the masses of the octet baryons as well as the electromagnetic and strong 
form factors of the nucleon. For different dressing functions we in particular 
compare the predictions for kaon photoproduction, jp — ► KA, and associated 
strangeness production, pp — > pKA, with experimental data. This leads to 
conclusions on the permissibility of different dressing functions. 
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I. INTRODUCTION 



The complexity of Quantum Chromodynamics (QCD) inhibits the computation of 
hadronic properties and reactions from first principles. As a consequence models that poten- 
tially imitate the essentials of the QCD dynamics have been developed in the past to describe 
hadrons. A relativistic description of baryons as three-quark bound states is provided by 
the solutions of the Bethe-Salpeter equations^ for quarks and diquarks which interact via 
quark exchange Once the full three-quark problem has been reduced to an effective 

two-body problem, the only model ingredients are the (di) quark propagators along with the 
quark-diquark vertices. It is hoped for that further progress in the study of the QCD quark 
propagator and two-quark correlations will eventually justify the reduction to quarks and 
diquarks in this approach to describe baryons. 

Actual calculations utilize either simplifying assumptions or phenomenological param- 
eterizations of the respective propagators and interaction vertices of quarks and diquarks. 
By choosing the simplest ansatze, i.e. free spin-1/2 and spin-O/spin-1 propagators for 
quarks and diquarks, respectively, various spacelike nucleon form factors have been success- 
fully reproduced 0. However, the naive use of perturbative (di)quark propagators leads 
to asymptotic states in the spectrum that carry the respective quantum numbers. Hence 
baryons would decay into quarks unless kinematically bound. This decay process would 
contradict the confinement phenomenon. In this paper we will therefore investigate the pos- 
sibility of incorporating confinement into the diquark model by suitable modifications of the 
quark and diquark propagators. Essentially these propagators are modified by multiplicative 
dressing functions to completely remove the poles that occurred in the perturbative propa- 
gators at the (di) quark masses. This enables us to calculate the spectrum not only of octet 
but also decuplet baryons f7j. Together with an ansatz for the quark-diquark bound state 
wave-function of the nucleon (Faddeev amplitude) such pole-free propagators have already 
been used to calculate nucleon form factors in the spacelike regime |||9|]. Unfortunately in 
this context the computation of the electro-weak form factors is not as simple as merely 
modifying the propagators. Since gauge invariance relates off-shell propagators and vertices 
it is obvious that any change in the propagators requires modifications of the effective inter- 



action with the electro- weak gauge bosons [10]. When incorporating gauge invariance in the 



model with free propagators the nucleon isovector magnetic moments come out too small 
by about 30% unless model parameters are used that do not properly reproduce the baryon 
spectrum ||. However, these unacceptable parameters result from requiring the A-isobar to 
be kinematically bound against its decay into free quarks. It is hoped for that when modeling 
confinement the results on the magnetic moments will also improve due to the modifications 
of the photon vertices which are mandatory when employing dressed (di) quark propaga- 
tors. A very serious disadvantage of the lack of confinement is that hadronic reactions with 
timelike momenta of the order of lGeV transferred to the nucleon, e.g. meson production 
processes, cannot be described properly. Again, the free-particle poles of quark and di- 
quark cause unphysical thresholds in these processes that are triggered by the poles in the 



1 For further details on the application of the Bethe-Salpeter formalism to QCD we refer to 
reviews [IT101 and references therein. 
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propagators. An appropriate modification of these propagators would not only remove the 
unphysical thresholds but also serve as an effective description of the strong interaction. 
Certainly, a relativistic description of such processes would be very desirable. At present, 
the covariant diquark model appears to be the only relativistic one that is both, feasible and 
applicable at this energy scale. 

As already mentioned we wish to eliminate the singularities associated with real timelike 
momenta in the (di) quark propagators that would lead to imaginary parts in those S-matrix 
elements that are calculated from diagrams containing internal quark loops. So, either these 
singularities are absent or their contributions cancel in some manner ||. The qualitative 
behavior described can be encoded in the following models (which are certainly not the only 
possibilities) for the quark propagator in Euclidean space, 



with 



S W ip) = %-^hK) , * = 0,...,3, (1) 
p z + mr q \rriq 



fo(x) = 1 (bare propagator) , (2) 

, , v _ 1 f x + 1 x + 1 \ 

h{X) ~ 2\x + l-t/d + x + l + i/dj ' [6) 

f 2 (x) = l-exp[-d(l + x)} , (4) 

f 3 (x, x*) = tanh [d (1 + x) (1 + x*)] . (5) 



The propagator ([3]) possesses complex conjugate poles JO]] such that corresponding virtual 
excitations cancel each other in physical amplitudes. Here m represents a parameter that 
would be interpreted as the quark mass if and only if the poles were on the real axis. 
In another scenario (f|), the dressing functions are chosen such that the propagators are 
entire functions and non-trivial in the whole complex plane [|l|J§. If they are required 



to be analytic, they must possess an essential singularity, at least for infinite arguments. 
Third, it might be helpful to approximate propagators by non-analytic functions @ and 
constrain them such that they asymptotically behave like l/|p| 2 for both, large spacelike and 
timelike momenta. Since we enforce the propagators to be free of poles, they must be non- 
analytic functions depending on both the particle momentum p and its complex conjugate p*. 
Consequently the quark-photon and quark-meson vertices are non-analytic and, even worse, 
translation invariance is lost in the solutions to the nucleon Bethe-Salpeter equation. These 
issues will be detailed in section [TI|. The trivial dressing function f corresponds to the bare 
propagator. We will consider this case only for comparison. 

In figure |l] we show fk(x) = fk(x)/(x + 1) for k = 0, . . . , 3 for real x. Note that these 
dressing functions are real in that case. We observe that fi(x) and fs(x) change sign (as in 
the case of a tree-level propagator) while the function f2{x) increases drastically. For asymp- 
totically large spacelike momenta the three model propagators S^ k \ (k = 1,2,3) match up 
with the bare propagator S^°\ Our present investigation focuses on the phenomenological 
implications of the so-modified propagators. 

This paper is organized as follows: In Section [TT] the covariant diquark model for baryons 
is presented. The corresponding Bethe-Salpeter equation that describes baryons as bound 
states of quarks and diquarks is derived in appendix [A[ The formalism of refs. Jl3||| for 
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FIG. 1. The propagator functions, fk(x) = fk(x)/(x + 1) for real x and for k = 0, ...,3, cf. 

eqs (|])-(||)- The thick solid line corresponds to the free propagator. Here we have set d = 1. 



calculating form factors is described for later determination of model parameters. Using the 
above given scenarios for implementing the confinement phenomena at the level of propaga- 
tors this will set the stage for the main topic of our paper: The sensitivity of the predicted 
observables on the various effective parameterizations of confinement. These parameteriza- 
tions concern the structure of the (di) quark propagators for complex momenta. In section [Tj I 
we will discuss the regime of complex momenta that is relevant for studying the baryon spec- 
trum as well as several production processes. In section IV we will describe the formalism 
necessary to compute various production processes in the diquark-quark model. These com- 
prise especially the cross sections for kaon photoproduction py — > KA and the associated 
strangeness production in pp — > pKA. We will proceed by presenting our numerical results 
in section |V|, including the determination of the model parameters. In particular we will 
compare the predictions that originate from the different dressing functions for the propa- 
gators. Finally, we will conclude by formulating criteria for phenomenologically permissible 
parameterizations of the propagators. Some derivations and technical details are relegated 
to four appendices. 



II. THE CO VARIANT DIQUARK MODEL FOR BARYONS 

A. The diquark quark Bethe Salpeter equation 

We start from the Faddeev formalism for three quarks and impose two essential as- 
sumptions to arrive at a Bethe-Salpeter equation that describes baryons as bound states of 
quarks and diquarks interacting via quark exchange. These assumptions are (i) all three- 
particle irreducible graphs can be safely omitted and (ii) the two-quark correlations can be 
approximated by separable correlations, the so-called diquarks. The actual derivation of 
the Bethe-Salpeter equation for the effective baryon-diquark-quark vertex functions (f) a is 
presented in appendix [A[ Using the definitions for total and relative momentum given as in 
figure ^ this integral equation reads 
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<t>tM P) = J2 [ 7S1 K tuM>>- I P) G% V {1 P) tfjril P) + (J — k) . (6) 
bb' J ^ 71 > 

Here Gj describes the disconnected quark-diquark propagator 

G%(1 P) = Sf V + I) D W {{1 - n)P - I) = Sf{l)D w {P - I) (7) 

Furthermore, the quark-diquark interaction kernel K contains besides the propagator of the 
exchanged quark also the diquark amplitudes xt defined via the separability assumption 

ti(kj,k k ;Pj,Pk) = ^2xi(kj,k k ) D a>al (kj + k k ) x?(Pj,Pk) (8) 

a, a' 

of the quark-quark t-matrix (cf. App. [AJ). The kernel explicitly reads 

Ktl P M I p) = xt Pl {i + vP, q) sr'(q) x b j , 7 'a(i,P+ vP) 

= xl^l, q )Sr\ q ) X l ya (q,p), (9) 

with 

q=(l- 2r])P -p-l = P-p-l, 

since P = Pi +pj +p k and 1 = 1 — rjP. The above relations also indicate the independence of 
the momentum partition parameter rj since the Jacobian of the transformation I —>■ I equals 
unity for fixed total momentum P. 

For the solution of the Bethe-Salpeter equation (||) we still have to choose the appropriate 
quantum numbers associated with baryons. This will be discussed in subsection [B] and we 
will find that the quark exchange (parameterized by the kernel K ab ) generates sufficient 
attraction to bind quarks and diquarks to baryons. For identical quarks antisymmetrization 
is required when projecting onto baryon quantum numbers. Fortunately, this does not alter 
the algebraic form of the Bethe-Salpeter equation (H). Rather, it simply implies that we 
may omit the single particle indices i on the quark propagators Si. Only when caring about 
the discrete quantum numbers we have to revert to these indices since they specify the 
summation order over color, flavor and Dirac indices in eq. (B|). Furthermore the functional 
forms of the diquark propagators D aa i and the vertices xl do not depend on the quark labels. 
These independencies are already indicated in eqs (0) and @ as we have omitted the quark 
labels for the momenta. 

In a self-consistent approach one would calculate the t— matrix from its own Bethe- 
Salpeter equation (A12). However, this is beyond the scope of the present investigation. 




FIG. 2. The coupled set of Bethe-Salpeter equations for the effective vertex functions <p a . 
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Instead we model the t-matrix by diquark correlators which have an analytic structure such 
that no particle interpretation for the diquark exists. We will restrict ourselves to the scalar 
and axialvector channels as these comprise the minimal set to describe octet and decuplet 
baryons. Furthermore these channels are generally assumed to be the most important ones, 
see refs. |1],|2| and references therein. The corresponding separable ansatz for the two-quark 
t-matrix reads 

t a w(h,h; Pl ,p 2 ) = x%(hP2) D(P 2 ) xUp,P 2 ) + X^(fc,P 2 ) D^(P 2 ) X^(p,P 2 ) • (10) 

Here we rewrite the diquark-quark vertices x 5 ^ as functions of relative,!; = ak\ — (l — a)k 2 , 
and total, P 2 = k\ + k 2 = pi + p 2 , momenta instead of the single quark momenta. In actual 
calculations, we choose for simplicity the symmetric momentum partition, i.e. a = 1/2. 
Shifting the value of a is possible, however, this complicates slightly the parameterization 
of diquark correlations, see the discussion below eq. ( |TTD and in ref. ||TB| . 



The diquark propagators in the scalar and the axialvector channel are modeled as 

D(P) = f ( , (11) 



P 2 + m 2 c \m 2 c/ 

^H-^^ + U-O^)/^). (12) 

The dressing function f(P 2 /m 2 ) is hereby chosen identical to the one for the quark prop- 
agator, i.e. either one of the forms (0,[|,|5]) . Note that the choice f(P 2 /m 2 ) = 1 and £ = 
corresponds to the free propagators of spin-0 and spin-1 particles. As a major purpose of 
the present paper we will study various deviations from the free propagators as an avenue 
to mimic confinement. In general, the dressing functions / are different in the scalar and 
axialvector channels as well as they are distinct from the one for the quark propagator. For 
simplicity, however, we will assume identical functions for all propagators. As we will not 
consider any axialvector diquark loops it is sufficient for the present purpose to use £ = 1, 
see ref. |7j where it has been shown that choosing £ = 1 leads to almost identical results for 
baryon amplitudes as £ = 0. 

If diquark poles existed in the t matrix, the diquark-quark vertices x an d X M would on- 
shell (P 2 = —m 2 c [ ax ]) correspond to diquark Bethe-Salpeter vertex functions. These vertex 
functions have a finite extension in momentum space and fall off fast enough to render all 
integrals finite. Empirically we assume that the corresponding scale is linked to the (inverse) 
proton radius. The conjugate vertex functions x are obtained by charge conjugation, 

x 5 ( Pl P) = C ( X 5 (- P ,-P)) T C T , (13) 
X?(p, P) = -C ( X »(-P, -P)f C T , (14) 

where T denotes the transpose. 

Let us now explicitly construct the vertex functions. They must be antisymmetric under 
the interchange of the two quarks. This entails 

xl [ % ] (P,P) = -xZ ] (-P>P) , , ■ (15) 
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Any two quarks within a baryon belong to the color antitriplet representation. Thus the 
diquark-quark vertices are proportional to the antisymmetric tensor €abd- Here A and B 
are the color indices of the quarks whereas D labels the color of the diquark. Furthermore 
the scalar diquark is antisymmetric while the axialvector diquark is symmetric in flavor. We 
maintain only the dominant components with regard to the structure in Dirac space.0 These 
are the antisymmetric matrix (7 5 C) for the scalar diquark and the symmetric matrices (7 M C) 
for the axialvector diquark. Considering, for the time being, only two flavors the vertices 
then read0 

4u(p, ^)U/2 = xlpiP) = 9s(l 5 CU V(p 2 ) ^ , (16) 

X^ip, P)U /2 = XUP) = 9a(rCU V(f) ^ . (17) 

Choosing the scalar function V to depend only on the squared relative momentum p 2 , these 
vertices are indeed antisymmetric with respect to exchange of quark labels for the partition 
a = 1/2. Otherwise a parametrization of V would depend on both p 2 and p ■ P in order 
to comply with antisymmetrization |13|J48| . However, complete independence for observable 



quantities on a could only be obtained by solving the Bethe-Salpeter equation ( |A12| ) for 
the two-quark t-matrix in which case the scalar functions V could depend on the quantity 
(p ■ P) 2 (for cr = 1/2) which is symmetric under quark exchange. In the actual calculations 
we will use a multipole form type ansatz 

v(x)=Vn{x)= ^y . (18) 



The overall strength of the diquark correlations given in eqs (p76|JT~7|) is governed by the 
"diquark-quark coupling constants" g s and g a . They could be determined by either imposing 
the canonical Bethe-Salpeter norm condition on x 5 ^ or by the solution to the differential 
Ward identity for the diquark-photon vertex which is sensitive to the substructure of the 
diquarks pj. For simplicity, we will fix g s from fitting the nucleon mass. When including 
axialvector diquarks we will assume the ratio g a /g s = 0.2 as suggested by the results of 
ref. || . In this manner the baryon Bethe-Salpeter equation (||) becomes an eigenvalue 
problem for the coupling constants g s and g a . 

Note that by parameterizing the quark-quark t-matrix we do not make any reference to 
the nature of the relevant quark-quark interaction. For example, to quantitatively include 
pionic effects we would have to solve Dyson-Schwinger equations for the quark propaga- 
tor and the Bethe-Salpeter equation with explicit pion degrees of freedom. Studies within 



2 The complete Dirac structure for the scalar diquark containing four independent tensors can be 
obtained by analogy from the one for pseudoscalar mesons. The complete Dirac structure for the 
axialvector diquark consists of twelve independent terms, four longitudinal and eight transverse 
ones. 

3 In the compact notation the indices a and of Xa/3 n °t only contain the Dirac labels but also 
those for flavor and color. 
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the Nambu-Jona-Lasinio model using diquark-quark correlations either in a soliton back- 



ground [IJJ or with explicit pion interaction between the quarks lead to a substantial 
gain in the binding energy. Since we determine the coupling constant g s from the nucleon 
mass, those studies suggest that the inclusion of pion degrees of freedom would merely lead 
to a shift of this constant. 

Equipped with the separable form of the two-quark correlations, see eq QlOD, and the 
functional form of the scalar and axialvector diquark correlations in eqs we will set 

up the effective Bethe-Salpeter equation for the nucleon. 

Upon attaching quark and diquark legs to the vertex functions <fi one obtains the Bethe- 
Salpeter wave functions ip. Equation ([]) can then be rewritten as a system of equations 
for wave- and vertex functions as defined in Appendix ||. Using the notations (|BTH B3|) we 
obtain 

Here G~ 1 (p, k, P) is the inverse of the quark-diquark four-point function which results from 
the quark exchange^. It is the sum of the disconnected part and the interaction kernel which 
contains the quark exchange, 

G-\p, k, P) = (2vr) 4 5\p - k) S-\ Pq ) o D-\p d ) (20) 

_\( -Xhvl) S T (g) xhvi) X^'ipi) S T ( q ) xV) 

2 V xhvl) S T (g) X^pD X"'(p2) S T ( q ) x"(p?) 

The flavor and color factors have been worked out and therefore x 5 {p 2 ) an d X M G° 2 ) from 
now on only represent the Dirac structures of the diquark-quark vertices (multiplied by 
the invariant function V n (p 2 ), cf. eq ( pj|) ). The freedom to partition the total momentum 
between quark and diquark introduces the parameter rj G [0, 1] with p q = i]P + p and 
Pd = (1 — 7])P — p. The momentum of the exchanged quark is then given by q = —p — 
k + (1 — 2i])P. The relative momenta of the quarks in the diquark vertices \ and \ are 
p 2 = p + k/2 — (1 — 3r])P/2 and p\ = p/2 + k — (1 — 3rj)P/2, respectively. Invariance under 
(four dimensional) translations implies that for every solution \l/(p, P]T}\) of the Bethe- 
Salpeter equation a family of solutions exists that have the form ^>(p + (772 — f]i)P, P; 772)- 
Considering the Bethe-Salpeter equation as a linear eigenvalue problem for \l/ (or $) in 
the coupling constant g s , translation invariance requires the coupling constant eigenvalue 
to be independent of r\ once a bound-state mass — P 2 = M 2 is fixed. This independence 
is exactly what one observes in the numerical solutions of the BSE, provided the analytic 
form of the dressing functions, eq ©— (P, is used. However, the ^-independence is lost when 
substituting non-analytic propagators such as those parameterized by the function / 3 , eq (0). 
Essentially the reason is that Cauchy's theorem does not apply to non-analytic functions. 
The difference in the eigenvalues of the Bethe-Salpeter equation under the variation of rj 
can be shown to equal a contour integral in the complex p-plane. This integral vanishes 
only if the integrand is an analytic function. However, when choosing d > 5 in eq (fj|), 



4 For convenience we have omitted the discrete labels. 
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the propagator resembles the free propagator in a large domain thereby mitigating the r\- 
dependence. 

The structure of the equations for the octet baryons is similar to that of the nucleon ([05). 



However, the number of Dirac structures $ 5 and $ M increases due to the possible different 
quark-diquark flavor configurations. These equations are given in full detail in refs. |7|,|T8 



Allowing for flavor symmetry breaking, that is induced by a difference between the masses of 
strange quark and up/down quark, discriminates vertex functions $ 5 and $ M with different 
diquark configurations j7|. As the A hyperon presently is of special interest, we list its three 
different correlations, 

$a~$! i; and $£. (21) 

Here, F\ = {d[us} — u[ds]} / y/2, F 2 = s[ud] and L = [d{us} — u{ds}}/\/2 refer to differ- 
ent quark-diquark flavor states. Antisymmetrized scalar diquarks are denoted by square 
brackets [. . .] and symmetrized axialvector diquarks by curly brackets {...}. Note that bro- 
ken S'f/(3)-flavor symmetry induces a component of the total antisymmetric flavor singlet 
4= [[stt]cf-|- [ud]s + [ds]u] into wave and vertex functions. In non-relativistic quark models 
with SU(6) symmetry such a component is forbidden by the Pauli principle, however, having 
non- vanishing lower components in the baryon bi-spinors does actually lead to such flavor 
singlet components. In actual calculations they turn out to be small |]7fl. 



B. Electromagnetic form factors 



To further constrain the model parameters, we calculate the electromagnetic form factors 
of the nucleon. In this section we provide the formalism and the corresponding results will be 
given in subsection V A . These form factors parameterize the nucleon matrix element of the 
current operator that describes the coupling of the photon to quark and diquark within the 
nucleon. Gauge invariance and the proper normalization of the nucleon charges are ensured 
when the current operator comprises all possible couplings to the inverse quark-diquark four- 
point function G~ x of eq ([Lj]) |l3|,P,[T8|j . The current operator is sandwiched between the 



Bethe-Salpeter wave-functions \l/ and \1/ of the final and initial state, respectively, according 
to Mandelstam's recipe p9| . In total, we have contributions from the impulse approximation 
which are described by the upper diagrams in figure ^| and contributions from the Bethe- 
Salpeter kernel which are given by diagrams of the type given in the lower part of figure |3|. 

To calculate the form factor diagrams, we need properly normalized Bethe-Salpeter 
wave-functions. This normalization is obtained from 



M N A 



d 4 p 

(2tt) 4 J (2tt) 



(,lk -Mk,p n ) 



P=Pn 



(22) 



where is the nucleonic bound state mass. The conjugated wave-function ^/ is in analogy 
with eqs ( |3~3|) and (|TJ) given by 



y(k,P n ) = r ld CV(-k,-P n ) T C 7 



(23) 



with rjd 
type. 



1 and rja = — 1 when the involved diquark is respectively of scalar or axialvector 
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Furthermore we need expressions for the photon vertices that appear in the diagrams. 
In ref. the seagull vertices describing the photon coupling to the diquark-quark vertices 
X have been derived for the scalar diquark. The coupling of the axialvector diquark to the 
photon has been studied in ref. |J. The photon vertices with quark and diquark must fulfill 
the differential Ward identities for zero momentum transfer to the photon 

r£(p^ = A ^(p,) , (24) 

OPq 

f3(p*p<) = ( r c? r°-) = W d D ~ 1(Pd) ' (25) 

Here comprises both photon vertices with scalar and axialvector diquarks. For conve- 
nience the discrete labels have been omitted. The Ball-Chiu construction of the longitudinal 
part of the vertices ensures that they obey both the differential Ward as well as the Ward- 
Takahashi identity. The latter identities constrain the vertices for finite Q. For the dressed 
quark propagators of eq (|IJ) the corresponding vertices finally read 



10 



rji(**,ft) = ~-f [l/fi(k 2 Jm q ) + l/Mp 2 q /m q )] 

i 1/ fAk 2 Jm q ) - l/Upl/m Q ) 
-o(^+P,r 9 5 2 - K^ + ^)-2^1 • ( 26 ) 

Z K q~ Pq 

This construction is valid for the analytic dressing functions /o, /i and /2. 

In the case of the non-analytic dressing function /3, we must specify the derivatives in 
eqs (P2]-p5|) with respect to the total bound state momentum P and the quark and diquark 



momenta p q and pa, respectively. We calculate form factors in the Breit frame, i.e. the 
temporal component of the momentum transfer is zero. Consequently, the relative momenta 
between the initial quark, k and the final diquark, p are real. They must be integrated over 
in the norm integral (p2|) and in the calculation of the diagrams of figure [| Let us consider 
the quark momenta which are defined as before: k q = r]Pf + k and p q = rjPi + p. We define 
the derivatives in equations fl22"D , (PH) and (E25|) as follows, 



d d d d , 

di^ = r? ^]' dp~ q 1 ^} = dp^\- [ } 

Of course, these are trivial identities when applied onto analytical functions. Derivatives with 
respect to the diquark momenta are defined accordingly. The nucleon charges obtained as the 
form factors at zero momentum transfer are then properly normalized. The corresponding 
proof utilizes the methods outlined in ref. [|TB|]. 

To comply with the Ward-Takahashi identity, the quark-photon vertex has to be modi- 
fied 

r£ 3 = ~Y (l/fs(k 2 q /ml,kf/m 2 q ) + l/fs(p 2 q /mlpf/ml)) 

- {k q + Pq r q q q t i [(h + h) - ^rn q ] 



2™ k\- Vq 



This vertex now depends on the four variables k q , k*,p q and p* and is also non-analytic as 
is the corresponding quark propagator S^ 3 \ The photon Ball-Chiu vertices with scalar and 
axialvector diquarks have to be modified using an analogous description ||18| . The coupling 
of the photon to the anomalous magnetic moment of the axialvector diquark and the vertex 
for photon-induced anomalous scalar-axialvector diquark transitions are transversal and 
need not be modified M. 



C. Strong Form Factors 

Here we consider the strong form factors g n NN and Qknk- These quantities are not only 
interesting in themselves but also enter the calculation of production processes like pj — > AK 
or associated strangeness production pp — > pKA. 

In figure |] we show the dominant contributions to the strong form factors. Here the meson 
directly couples to one of the baryon constituents. Keeping only such direct couplings while 
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FIG. 4. Dominant diagrams for the strong form factors g w NN an d 9k N A- The incoming proton 
carries the momentum Pj while Pf is associated with the outgoing baryon. The incoming meson carries 
the momentum Q and couples to the quark (left panel) or to the diquark (right panel). 



omitting those to the exchanged quark defines the impulse approximation that we will adopt 
hereQ. The two diagrams shown in figure ^ actually correspond to a number of diagrams 
which differ by the type of the involved diquarks. Let us first consider the process in which 
the meson couples to the quark. For Qknk only one possibility exists: the diquark has to be 
a scalar wd-diquark since this is the only overlap between the wave-function of the proton 
and the wave-function of the A. For g^NN both scalar and axialvector diquarks need to be 
taken into account. The second important contribution represents the coupling of the meson 
to the diquark. For the diquark part we do not have to distinguish between g^NN an d 9kna- 
That is, in both cases the diquark associated with the momenta p + or p_ may be scalar or 
axialvector. 

The meson-quark vertex is the solution of a separate Bethe-Salpeter equation which 
has been extensively studied, see and references therein. In the chiral limit this Bethe- 
Salpeter equation becomes formally identical to the Dyson-Schwinger equation for the scalar 
self energy function B(p 2 ) when only the leading Dirac structure is considered, i.e. 



Q 



; ,_ , -r m (P-,P + ) = ^jl 5 {B(pl) + B(pl)} . (29) 

Here / is the meson decay constant. 

The structure of the meson-diquark vertices is constrained not only by Lorentz covari- 
ance and parity but also by the Bose-statistics for the two involved diquarks. We thus 
parameterize the pseudoscalar meson axialvector diquark vertex as 

r£=-i^je pX » u (p-+P + rQ». (30) 

Here the superscripts p, A denote the Lorentz indices of the incoming and outgoing ax- 
ialvector diquark, respectively. The nucleon mass M has been introduced to define the 



5 Contributions to the nucleon electromagnetic form factors beyond the impulse approximation 
that arise from the coupling to the exchanged quark have been thoroughly discussed in refs [| , 13,18 1 . 
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dimensionless coupling constant n aa . Furthermore m is the average of the masses of the 
constituent quarks in the diquarks. The corresponding ansatz for the scalar-axialvector 
transition reads, 

77? 

K a = -K sa -Q<>, (31) 

where the definitions are those of eq ( j3Tj) and n sa is again a dimensionless constant specifying 
the overall strength of the vertex. The vertex fl3~l|) describes the coupling of the diquarks 
to the derivative of the pseudoscalar mesons. Such a construction is suggested by the chiral 
structure of the strong interactions that can be written as expansion in the derivatives of 
the Goldstone bosons, at least in the chiral limit. 

Having collected all ingredients we may now proceed and compute the diagrams in fig- 
ure |j. According to the Mandelstam formalism jn| the diagram shown in the left panel 
translates into an expression of the form 

/ -0jI^ ) f( ( lf^f) S (P+) T rn(p-,P + )S(p^p(q,P l )D(Pci), (32) 

where we only indicated the general structure, i.e. we omit indices that are associated to 
the coupling and propagation of axialvector diquarks. The conjugated vertex function, <§ 
relates to the vertex function, $ as the conjugated wave function (p3|) to the wave function: 

$(p, P) = VdC $(-p, -P) T C T (33) 

with T]d = 1 and rjd = — 1 when the involved diquark is respectively of scalar or axialvector 
type. We note that $(p, P) also solves the Bethe-Salpeter equation (|). We denote the loop 
momentum by q and introduce the momentum partition, 

p_ = q + rjPi, p + = p- + Q = q f + rjPf and p d = -q + (1 - T))Pi . (34) 

Again, r) G [0, 1] is the momentum partition parameter. For the diagram in the right panel 
quark and diquark propagators need to be exchanged. 



III. BOUND STATE REACTIONS AND KINEMATIC AL CONDITIONS FOR 

COMPLEX MOMENTA 

In this section we will discuss that regime in the complex momentum plane where we 
need to know the quark and diquark propagators in order to solve the Bethe-Salpeter equa- 
tion (|i~9|) and compare that regime to the one that enters the computation of the production 
processes like pry — > AK and pp — > pKA. In principle these propagators can be calculated 
using Dyson-Schwinger equations and also respective lattice measurements should be 
available in the near future; for preliminary results see e.g. refs. [^,^T[|. Both methods 
comprise the non-perturbative dynamics and should therefore give the basic ingredients to 
describe hadrons as bound state of quarks. However, both approaches are set up in Eu- 
clidean space and one has to revert to extrapolations when the propagators are demanded 
for timelike momenta. If we wanted to perform an appropriate analytic continuation from 
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FIG. 5. The complex g 2 -plane. The interior of the parabola is needed for the calculation of the 
diagram in the left panel of figure ^, x is defined as x = (rjM + E) 2 . Note that in the case of the 
Bethe-Salpeter equation we substitute x = i] 2 M 2 . 



Euclidean back to Minkowski space we would even require the propagators in a region of the 
complex momentum plane]]. In order to calculate amplitudes of physical processes between 
on-shell particles using the Euclidean Bethe-Salpeter formalism the temporal components 
of the external momenta must be purely imaginary. In this framework the momenta become 
complex. In these calculations therefore the structure of the propagators in the complex 
momentum plane is essential. Furthermore, it is important for the phenomenological pa- 
rameterization of confinement. As repeatedly mentioned we comprehend the confinement 
phenomenon as the absence of poles on the timelike g 2 -axis in the propagator of colored 
"particles" . 

The Bethe-Salpeter equation is most conveniently solved in the rest frame of the bound 
state, P = (0,iM). Here we want to specifically discuss the kinematical domain that is 
probed by the (di)quark propagators in the bound state rest frame. In eq (]TQ|) the loop 
momentum, k, relative between quark and diquark is chosen to be real. Hence the temporal 
component of the quark momentum k q = rjP + k becomes complex. The values of k 2 that 
are covered when integrating over k lie within a parabola that opens towards the spacelike 
axis, cf. figure || The intercept of the parabola with the real axis is at (small) timelike 
kg = —(t]M) 2 . Thus, solving the Bethe-Salpeter equation mainly probes the behavior 
of the quark propagator for spacelike momenta. Since mainly the spacelike momenta are 
relevant, the propagators that are parameterized by the dressing function (|4]) approach the 
free propagators in the limit d — > +00. In this limit the dressing functions (H) and (|S|) 
approach the bare propagators, both in the spacelike and timelike regions. 

Next we will explore the q 2 regime relevant for production processes like kaon photopro- 



6 In the Dyson-Schwinger approach the corresponding integral equation should be used for this 
analytic continuation. Relying on a numerical solution that is only known for a finite set of 
Euclidean momenta is not sufficient because its analytic continuation away from that set cannot 
be determined. 
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FIG. 6. Left panel: Main contribution to kaon photoproduction p'j — > KA, right panel: Handbag 
diagram contributing to the reaction pp — ► pKA as a subprocess. The incoming pion couples to the 
'spectating' proton. 

duction. The contribution to the reaction jry — > AK that involves a quark loop is shown 
in the left panel of figure §. It turns out that it suffices to consider a parabola shaped 
region of the complex g 2 -plane (i.e. it is sufficient to consider only this momentum and 
ignore the others). This can be understood in at least two ways: We could use the wave- 
functions rather than the vertex-functions for the calculation of the diagram. In this case 
the propagators that depend on p q , k q and pd are included in the wave-functions and there 
would be no necessity to treat them separately. Nevertheless, considering the propagators 
S(p q ), S(q), S(kq) and D(p d ) separately we find that among all the internal momenta in the 
diagram it is q that reaches farthest in the timelike regime. Thus the following analysis for 
p q and k q would yield less restrictive conditions. 

From the momentum routing shown in the left panel of figure ^| we have: 

q = P- + p 7 = rjP + I + p y (35) 

where i] is the momentum partition parameter (p_ = 77 P + I, pa = (1 — rj)P — I) while / 
refers to the loop momentum. We chose the loop momentum to be real which implies that 
the external momenta like P and p 7 must have an imaginary temporal component in order 
to correspond to physical particles. For the following kinematical considerations we choose 
the proton rest frame and take the photon to propagate along the y-axis, 

P = (0, %M) , p 7 = (0, E, 0, iE), 1 = {I, h ) (36) 

Hence the momentum entering the quark propagator becomes 

q 2 = {-rfM 2 + l 2 - 2t]ME + 2El y ) + % (2r]M + 2E) l 4 , (37) 

where the real and imaginary parts of q 2 have been separated. This shows that we need 
to know the propagator S(q 2 ) at complex q 2 in order to be able to compute the handbag 
diagram shown in the left panel of figure || The set of values of q 2 that occur has already 
been shown in figure [5]. The situation seems to be completely parallel to what we found for 
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the Bethe-Salpeter equation; in both cases we need to know the propagators in a parabola 
shaped region of the complex plane. The intercept with the imaginary axis is in both 
cases minus two times the intercept with the real axis. However, there is one important 
difference. For the production processes the intercept with the real axis does depend on the 
photon energy E, more precisely: — x = —(rjM + E) 2 . Thus for E — the computation of 
the handbag diagram shown in figure |6] uses the same region of the complex plane that is 
necessary to solve the Bethe-Salpeter equation. However, for E > the parabola is shifted 
in the direction of the negative real axis. 

The threshold for kaon photoproduction is at E slightly less than lGeV and the cross 
section has been measured []22] up to E ~ 2GeV. This implies that the handbag diagram 
'probes' the quark propagator much farther into the timelike region than the Bethe-Salpeter 
equation. 

The second production process we are especially interested in is associated strangeness 
production, pp — > pKA. This reaction can be described similarly to the standard picture of 
the nucleon-nucleon interaction by one-boson exchange. That is, one of the incoming pro- 
tons acts as a meson source and the emitted off-shell meson couples to one of the constituents 
of the baryon; the corresponding subprocess is shown diagrammatically in figure || 

The analogous kinematical analysis for strangeness-production exhibits the same quali- 
tative features. That process as well 'probes' a parabola shaped subset of the complex plane, 
whereby the parabola is somewhat broader than the one in figure [5| However, there again 
is an important difference: the parabola does extend only up to q 2 « -0.53GeV 2 into the 
timelike region. That is, the reaction pp — > pKA 'probes' the propagators in essentially the 
same region as the Bethe-Salpeter equation does. It is therefore not as sensitive as kaon 
photoproduction to the behavior of the propagator in the timelike region. 

The main conclusion of the above discussion is that certain production processes may be 
significantly more sensitive to the structure of the (di) quark propagators than the Bethe- 
Salpeter equation and thus the baryon spectrum. Hence the study of such processes should 
provide important information about these propagators. 

IV. PRODUCTION PROCESSES IN THE DIQUARK QUARK MODEL 

In this section we present the key issues of the formalism to compute the cross sections for 
kaon photoproduction, 7p — > AK and the associated strangeness production, pp —>■ pKA. In 
the diquark-quark model relatively few diagrams contribute to these processes and therefore 
we may analyze these reactions in detail. For further details on the definition of the involved 
observables and the relevant kinematics we refer the reader to appendix |C[ 

As already indicated in the discussion of the strong form factors we consider the pseu- 
doscalar mesons as additional model degrees of freedom. This does not imply any double 
counting because the model interaction (diquark exchange) does not lead to bound (would- 
be) Goldstone bosons. Thus we also include intermediate pseudoscalar mesons at tree level 
when computing the above mentioned observables. The relevant diagrams are shown in fig- 
ures [7] and H As a general remark we note that these diagrams need to be computed in any 
covariant diquark-quark model. However, the propagators that are essential components 
of these diagrams are specific to our model, cf. eqs (0)-©. Furthermore, the covariant 
wave- or vertex-functions that also enter these diagrams are obtained as solutions of the 
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FIG. 7. Main contributions to kaon photoproduction jry — > AK. The incoming proton and the 
outgoing A carry the momenta P and Pa respectively. The lower part of the figure shows the tree level 
diagram that models the exchange of a virtual kaon. 

Bethe-Salpeter equation. Since this equation is subject to the model propagators they enter 
the calculation not only explicitly but also implicitly. 

A. Kaon Photoproduction 

In this subsection we will discuss kaon photoproduction py —> AK within our covariant 
diquark-quark model. Some more details and technicalities of the calculation are given in 
appendix |C 1| . 

We show the dominant diagrams in figure [7]. The internal momenta of the (uncrossed) 
'handbag diagram' are defined according to 

p_=p + T) p P, p d = -p + (1 - r) p )P , q = p„+p^, 

p+ = q-pK, PA = P + P-y~PK, Pf = P+ (1 - n P )P - (1 - Va) f a • (38) 

Here rj p and r]\ are the momentum partition parameters of the proton and the A, respectively. 
Both, rj p and r/\ can be chosen independently in the range < r) p , T]x < 1- 

The two 'handbag diagrams' model the coupling to one of the constituents. They are 
calculated within the Mandelstam formalism. This yields 
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diagram # 


1 


2 


3,5 


4 


diquark content 


s 


s,a 


cLS£l ? clclS, cLclcl^ SclS ; Sclcl 


SclS, Sclcl 



TABLE I. The diquark content of the diagrams shown in figure || with the numbers referring to 
the specific diagram. Here 's' and 'a' indicate scalar and axialvector diquark respectively. A sequence 
with three entries runs clockwise and starts at the incoming proton. 



Ai = if {$ A (p f ,P A )S(p + )r K (q,p + )S(q)} {T,(p_,q)S(p_)<S> P (p,P)D(p d )} (39) 

for the amplitude of the uncrossed handbag diagram. Here $a and $ p are respectively the 
vertex-functions of the A and the proton as discussed in section [TI]. Furthermore Tk is 
the meson-quark vertex that has been discussed in the preceding subsection. The photon- 
quark coupling, T 7 is described by the Ball-Chiu vertex |23[] or its generalization to the case 
of non-analytic propagators, see section [H| The Ball-Chiu vertex has been constructed to 
satisfy the Ward identity. It reduces to the bare vertex in the limit that both momenta p and 
q are large. The Ward identity constrains only the longitudinal part of the vertex. Various 
ansdtze for the transversal part of the vertex have been proposed (cf. ref [0 and references 
therein). While those ansdtze solve problems related to multiplicative renormalizability 
and gauge invariance, the transversal part is generally assumed to be of minor influence 
on the resulting cross sections. Thus we will henceforth neglect the transversal part of the 
quark-photon vertex. Although the form of this vertex is not model specific, it contains 
the self-energy functions and thus it implicitly depends on the model propagators. The 
expression for the crossed handbag diagram can be easily inferred. The tree level diagram 
models the exchange of a virtual kaon and is expected to yield a non-negligible contribution 
for large photon energies. For the photon-meson coupling we use a bare vertex multiplied 
with the kaon electromagnetic form factor (see Appendix |C 1| ) while the meson-baryon vertex 
is proportional to gKNA^Q 2 ) that has been discussed in subsection p C . 



The 'handbag diagrams' shown in figure [7| probe the propagators not only for space- 
like momenta but also for comparably large timelike momenta, as we have emphasized in 



section |TJ. This sensitivity to the behavior of the propagators for timelike momenta distin- 



guishes the reaction pr/ — > AK from most other production processes. 



B. Associated Strangeness Production 



Here we apply the covariant diquark-quark model to associated strangeness production, 
pp — > pKA. Again, some technicalities are relegated to appendix |C 2| . We describe the 
reaction pp — > pKA as a sum over one-boson exchange contributions for which we consider 
the exchange of pions and kaons. The main contributions are shown in figure |8|. Other 
diagrams like planar kaon exchange or crossed pion exchange with couplings to the quark 
are excluded by the flavor algebra. We neglect diagrams in which the exchanged particle 
and the emitted kaon couple to different constituents, because they imply a large relative 
momentum at the baryon vertex which is strongly suppressed. In addition we omit the 
direct coupling of the pion to the emitted kaon. 
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FIG. 8. Main contributions to associated strangeness production pp — > pKA. The momenta of the 
initial protons are given by Pi, i = 1,2. The momenta of the final proton and A are denoted by Pf and 
Pa, respectively, px refers to the momentum of the final kaon. The diquark content of the diagrams 
is listed in Table |. 



19 



Similar to the calculation of gxNA in section fTO, the flavor content of the A prohibits 



an axialvector diquark in the pion exchange diagram However, this does not apply 
to the kaon exchange diagram 7^2 as the diquark mediates between two protons. Due to 
parity conservation we only have to consider scalar-axialvector and axialvector-axialvector 
diquark transitions at any meson-diquark vertex. Flavor algebra also shows that the axi- 
alvector diquark components of the incoming proton cannot contribute in diagram #4. The 
diquark content of the diagrams shown in figure ^| is listed in table |. Counting the diquark 
combinations in table | we arrive at 15 diagrams. This number is actually doubled because 
all diagrams have to be antisymmetrized with respect to the two incoming protons. 

As an example, we outline the calculation of one of the two planar pion exchange dia- 
grams. The calculation of the other diagrams is very similar. The amplitude of diagram 
can be factorized according to 



m=c \wt*l) h ' (40) 

where £ denotes the form factor part and 7i denotes the loop part of the diagram. The 
mass and the momentum of the intermediate pion are denoted by m n and Q, respectively. 
The factor C essentially equals g w NN, 

£ s ,s> = u s t(P f ) i^g^NNiQ 2 ) u s (P 2 ), (41) 

with spinor indices s and s'. In the remaining 'handbag part' 7i the conventions for the 
loop momenta can be extracted from diagram #1 of figure [8]. Essentially they are given in 
eq ( |55D with the substitution p 7 — > Q and similarly for the baryons. The 'handbag part' of 
the pion exchange diagram can then be written as 

n = i J -0- {$a(p/,Pa) S(p + )T K (p + ,q) S(q)} {TM1P-) S(j>-)D{p A )* P faP x )} , (42) 
where isospin as well as Lorentz indices have been omitted for simplicity. 



V. NUMERICAL RESULTS 

After having outlined the model calculation we are now prepared to present the numerical 
results. Here we focus on studying the effects of the different model propagators (0)-(§) on 
the calculation of and the predictions for the above mentioned processes.[] As mentioned 
earlier, this is the main purpose of the present study. Numerical results for the form factors 
obtained with the tree level propagators (0) can be found in ref. 0. 



7 A few selected numerical results using other forms of dressed model propagators have been 



published in refs. |24|. 
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A. Parameter Fixing: Masses and Electromagnetic Form Factors 



We fix the model parameters, see table O, from the octet baryon masses and the nucleon 
magnetic moments. The numerical details for solving the octet baryon Bethe-Salpeter 
equations and the computation of the form factors are thoroughly discussed in ref. |18| . 

Within the required numerical accuracy we have assured the above described indepen- 
dence of the octet masses of the momentum partition parameter r\ when analytical propaga- 
tors are used. As argued before, this invariance does not hold for non-analytic propagators. 
In these cases we chose r\ to be close to its non-relativistic value m q j {m q + m d ) where m q 
and md denote quark and diquark mass parameters of the flavor channel associated with 
the considered baryon. This choice is natural since other ones yield larger eigenvalues of the 
Bethe-Salpeter equation. We take the physical nucleon mass to fix the scalar diquark cou- 
pling g s and the A mass to determine the strange quark mass parameter m s . By reproducing 
the phenomenological dipole fit for the proton electric form factor, Ge we essentially fix the 
diquark width A. Subsequently we are enabled to compute the proton and neutron magnetic 
moments, fi p and fi n as well as the masses of the £ and H baryons. For that calculation we 
assume isospin symmetry, m u = m d . 

In table |I| we list the six parameter sets that we will employ to compute the strong form 
factors and observables of production processes later on. The first four sets are restricted 
to the dominant scalar diquark correlations. In set I we consider the pole-free exponential 
dressing function, f 2 , while the sets II and III are associated with dressing functions of the 
Stingl type, f±. These two sets differ by the value of d that characterizes the separation of 
the complex conjugated poles. Finally set IV assumes the non-analytic pole-free dressing 
function, / 3 . As already indicated the dressing of the propagators increase the predicted 
proton magnetic moment when all other model parameters remain unchanged. Using the 
parameters of set II but free propagators yields fi p = 2.27 while the Stingl-type propagators 
result in fi p = 2.46 and \x v = 2.64 for d = 8.0 and d = 4.0, respectively. The magnetic 
moment of the proton falls a little short for the sets II and IV. The overall picture emerges 
that the restriction to only scalar diquark correlations produces too large ratios |/i n //z p | and 
rather large mass splittings between the octet baryons, especially between E and A. 

As shown in figure |9| all sets reasonably well reproduce the electric form factor, Ge- 
Our results are confined within a region that is characterized by less than approximately 
15% deviation from the dipole fit. This deviation, although rectifiable by refining the time- 
consuming parameter search, is of no significance for the conclusions that we will draw from 
our results for the production processes. This will become clear from the discussions in 



section VD 



The calculations with the parameters sets V and VI include a moderate admixture of 
axialvector diquarks, g a /g s = 0.2. For simplicity the axialvector diquark masses are chosen 
identical to the scalar ones. Here we particularly consider the Stingl form, f\ (set V) 
and the non-analytic form, / 3 (set VI) since later we will find that the exponential form, / 2 
produces unacceptable results for the production processes. Upon inclusion of the axialvector 
diquark the good description of Ge remains unchanged while the ratio |/i n //x p | and the mass 
splitting between £ and A even improve. For set VI the predicted octet masses are almost 
indistinguishable from their experimental values. As already observed in ref. || and as is 
exhibited in the right panel of figure 0, the ratio Ge/Gm calculated with axialvector diquarks 
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TABLE II. The six parameter sets of the model investigated here and the respective results for 
the nucleon magnetic moments and the octet masses. Calculations using the first four sets involve only 
scalar diquarks, whereas the sets V and VI also include axialvector diquarks. The parameter ( determines 
the diquark mass parameter (scalar and axialvector), nid = ((m a + mb), with the mass parameters m a ,b 



of its constituent quarks. The parameter A determines the width of the diquark amplitudes, see eq (18). 
For set I the corresponding shape of the amplitudes was chosen to be a quadrupole (n = 4), for the 
other sets we fixed it to be a dipole (n = 2). 



included comes considerably closer to the experimental values than in a calculation that 
omits these degrees of freedom (sets I-IV). As explained in ref. ||, increasing the strength 
of axialvector correlations in the proton forces the ratio Ge/Gm to bend to lower values. 
This also suggests that in order to precisely reproduce the empirical result we would need 
an even slightly larger axialvector coupling than the assumed g a /g s = 0.2. 

All sets predict the maximum of the neutron electric form factor to lie between 0.025 
and 0.04. This is only about half the value extracted from recent experiments ]27 , |28| . This 



form factor is a result of delicate cancellations between the contributions from the individual 
diagrams shown in figure [| Hence it is quite sensitive to small changes in the parameters. 
Within this model approach improved descriptions for this form factor can be found in 
refs. [IBI-pl]. 



In a previous study || that employed free quark and diquark propagators we were unable 
to reproduce the nucleon magnetic moments and the A mass simultaneously. The kinemati- 
cal binding of the A required a large constituent quark mass, m q = 0.43 GeV, which in turn 
decreased the magnetic moments (in magnitude). Furthermore the use of free propagators 
enforced moderate axialvector diquark contributions (about 25%) to properly describe the 
ratio Ge/Gm of electric and magnetic form factors for Q 2 up to 2GeV 2 . In contrast, the 



introduction of dressing functions for the quark-photon vertex (|26|,|28|) allows us to choose 



rather large up quark mass parameters around m u = 0.45 GeV and still obtain a proton 
magnetic moment that agrees with experiment reasonably well. 

Let us briefly reflect on the accuracy of our calculations. Due to the Monte-Carlo 
integration of the diagrams given in part b) of figure (with 7.5 ■ 10 5 grid points for sets I-IV 
and 4 ■ 10 5 grid points for sets V and VI) the absolute numerical error for \i v is 0.02 and for 
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FIG. 9. Left panel: Proton electric form factor normalized to the dipole fit. The experimental 
data are taken from ref. Right panel: The ratio {h p Ge)/Gm for the proton with the experimental 
data published in ref. [26|. 



fi n it is 0.03. The statistical error for the electric form factor is below 0.002 up to momentum 
transfers of 1.7GeV 2 . For the sets IV and VI a systematic relative error in the electromagnetic 
form factors is found that increases slowly to about 5 % at Q 2 = 1.7GeV 2 . Furthermore, 
we used an expansion in Chebyshev polynomials for wave and vertex functions because 
this expansion can be unambiguously continued to complex arguments. Upon employing 
analytical propagators we have obtained identical results when using the approaches wherein 
either the wave functions \l/ or the vertex functions $ must be continued analytically ||13 



For non-analytic propagators, the corresponding solutions of the Bethe-Salpeter equation 
are non-analytic as well. Thus the analytic continuation produces errors which can be 
estimated by comparing the calculations using wave functions or vertex functions. In order 
to restrict the resulting discrepancy to below 5%, we have to choose d > 5 and rather large 
quark and diquark mass parameters. This then leads to nucleon solutions of the Bethe- 
Salpeter equation that resemble an analytic function in the kinematic domain needed for 
the ongoing computation. 



B. Strong Form Factors 

Above we have utilized baryon properties to determine the model parameters that en- 
ter the Bethe-Salpeter problem. Furthermore the meson-quark vertex is governed by the 
appropriate Ward-Takahashi identity. Finally the meson-diquark coupling constants K sa 
and K aa in eqs (|30]) and (|31|) have already been determined in ref. ||. Thus we are now 



completely prepared to compute the loop-integrals like that in eq ( |3"2"D and that appear in 
figure |. Subsequently we may extract q^nn and gKNA- 

There have been numerous experimental efforts to determine the strong form factors. 
An extended discussion of the phenomenological value of g n NN and a comprehensive list of 
related references is provided in ref. |29[| . Mainly the quoted discrepancies are subject to 
different analysis of available data. For the purpose of the present work it is sufficient to 
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know that the quoted data are of the order g^NN ~ 14. Unfortunately the measurements 



of 9kna have not yet reached a satisfactory accuracy. The authors of |30[ have extracted 



\9kna{Q 2 — — = 13.7 ± 0.9 from the LEAR-data; but other analysis have partially 



yielded quite different results |31]j32 



In figure [H] we display the numerical results for the g n NN(Q 2 ) and Qkna(Q 2 )- Both form 
factors have been calculated in the Breit-frame. This frame is peculiar because for different 
masses of the initial and final baryons a numerically save treatment induces a lower bound 
(Q 2 > A| F = M\ — Mjf « 0.4GeV 2 ) for the momentum, Q, of the coupling meson. In the 
case of gxNA we have fitted the computed form factor to rational functions and extrapolated 
those functions to Q 2 — > 0. The resulting coupling constants are shown in table fTJ. For 
the special case Q 2 = we have verified that this treatment yields the same result as a 
calculation in the lab-frame. 

For both form factors we observe a qualitative difference between calculations with or 
without axialvector diquarks included. Whereas for all parameter sets with only scalar di- 
quarks the computed pion-nucleon form factor very well reproduces the experimental data, 
we find that for those sets that include axialvector diquarks the computed form factor over- 
estimates the data. This could be due to the omission of subdominant amplitudes in the 
meson-quark vertex. We performed a rough estimate of the influence of the first subleading 
amplitude by using a simple parametrization and indeed found negative corrections to the 
pion-nucleon form factor of about 30%. Future calculations should include these contribu- 
tions in a selfconsistent way. On the other hand in each of the subsets (I-IV) and (V-VI), the 
absolute value of the couplings at Q 2 = and the respective slope are almost independent of 
the parameter sets and even of the propagator type. To further analyze the structure of the 
form factors we have disentangled the various contributions in figure 11. This figure shows 
that at small positive Q 2 the contribution from the coupling of the meson to the scalar quark 
is clearly dominating, whereas for larger momenta the diquark contributions take over. 

As the form factors serve as input for later calculations we have conveniently fitted our 
numerical results to rational functions allowing, however, for a non-integer exponent, 

g*B»m = 9<t>BB' (A| F ) (^fr) P • ( 43 ) 

The result of this procedure are summarized in table [TTI|. To compare with results of other 
model calculations we have additionally fitted our results to monopole form factors. This 
has yielded scales A between 200 — 300MeV for the parameter sets with only scalar diquarks 
and scales around 500MeV for the sets including axialvector diquarks. As can be seen from 
figure |TT| the additional contributions fall substantially slower than the scalar diquark one 
and become dominant for large Q 2 . This effect can be interpreted as 'hardening' of the form 
factor. In agreement with the results from ref. || our form factors are much softer than 
those usually substituted in one-boson-exchange potential models for production processes. 
Those empirical scales for the monopole form are larger than 1300MeV f33fl . However, other 
theoretical approaches, e.g. lattice measurements or QCD-sum rule calculations indicate a 
monopole behavior with much smaller scales, cf. ref. ]34| and references therein. 

Our prediction for gxNA^Q 2 = 0) is comparable to those found of QCD-sum rule 
or Skyrme model calculations but somewhat smaller than the chiral bag model result, 
cf. ref. and references therein. Extrapolating our gxNA{Q 2 ) to the kaon mass shell 
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diquark: 
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9ttNN 

exponent 
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1 "3 A 
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14 ^ 
5.1 
1132 


i a n 
5.2 
1252 


18.2 
1.4 
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17.5 
1.6 
778 


1 ^ 4 


9KNA 

exponent 
scale [MeV] 


7.98 
10.6 
1786 


7.39 
5.9 
1368 


8.25 
5.8 
1391 


8.12 
6.4 
1554 


11.97 
1.3 

642 


10.23 
2.0 
950 





TABLE III. Numerical results for the absolute values of the couplings g n NN and gxNA at zero 



squared momentum and for the best rational fits to the curves. The parameter sets are defined in 
table |Tj. The entries 'exponent' and 'scale' refer to the variables p and A in the fit (ft3|). 



Pion Nucleon Form Factor Kaon Nucleon A Form Factor 




FIG. 10. The model prediction for the strong form factors g n NN and gxNA- The parameter sets 
(l-VI) are defined in table |j. For gKNA the results in the region Q 2 < A| F are obtained from a rational 
fit (see text). 



Q 2 = —M\ yields values in the range 16.3 < Qkna < 19.3. This is slightly above the 



ballpark of the numbers extracted from experiment [3D . 

The comparison between g n NN and gxNA suggests three different scenarios of SU(3)- 
flavor symmetry breaking that are illustrated in figure |T2|. The most obvious symmetry 
breaking stems from different quark mass parameters m u ^ m s . Secondly, due to the flavor 
algebra the process in which the axialvector diquark acts as a spectator, only contributes to 
g-nNN- Finally there are different decay constants fx ^ fn that factorize in the meson-quark 
vertices. These three effects cause gxNA(Q 2 ) < g-nNNiQ 2 ) independently from the type of 
the propagator or adopted model parameters. We see from figure |12] that at moderate and 
large Q 2 the different decay constants dominate the symmetry breaking effects. For small 
Q 2 the mass differences are essential. 
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Pion Nucleon Form Factor Kaon Nucleon A Form Factor 




FIG. 11. The distinct contributions to the form factors g n NN and gxNA calculated using parameter 
set V. 



Flavor Symmetry Breaking 

decay constant vs. strange quark mass (set VI) 

20 | , , , , , , . 




FIG. 12. Analysis of flavor symmetry breaking in gxNA- 
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FIG. 13. Total and differential cross section of the process pp 
section (right panel) is considered at an excess energy of 138MeV. 



pKA. The differential cross 



C. Associated Strangeness Production 

The COSY-TOF collaboration has measured the cross section and the polarization for 
pp — > pKA at 55MeV and 138MeV above threshold [[JIJ. In addition there are also data for 
the depolarization tensor D^n from the DISTO collaboration at SATURNE II at an excess 
energy of E 



431MeV [37 



This tensor is an especially interesting observable. Eventually 
Dnn might provide further information on the spin structure of the nucleon because it 
describes the transport of spin from the initial to the final states (cf. appendix C 2| for 
appropriate definitions). 

Our numerical results for associated strangeness production are shown in figure [13]. The 
comparison with the empirical data clearly shows that the propagator with an exponential 
dressing function (set I) yields unacceptable results. As discussed in section |T| the mecha- 
nism is that by increasing the beam momentum larger timelike momenta appear in the loop 
propagators and hence the cross sections suffer an exponential enhancement. This effect is 
most strongly pronounced at forward and backward angles in the differential cross section. 

All propagators that do not involve the exponential dressing function underestimate the 
cross section for pp — > pKA considerably. Only the parameter set VI can be considered 
to be at the right order of magnitude. Generally we find that the inclusion of axialvector 
diquarks improves the agreement with the data. This is not only the case for the total cross 
sections but also for the shape of the differential ones. The huge dip that arises for the 
parameter sets I- IV at directions perpendicular to the beam axis is considerably damped by 
the axialvector diquark contributions although it is still too deep. 

The two distinct contributions to the cross section that can be characterized as being 
associated with pion or kaon exchange (cf. figure |S|) lead to significant interference cancella- 
tions for the depolarization tensor D NN . The kaon exchange processes generate the outgoing 
A in the form factor part of the diagrams. Here the 75-structure of the kaon vertex leads to 
a spin flip from the incoming proton to the A because of parity conservation. Therefore kaon 
exchange diagrams provide negative contributions to the polarization tensor. In the pion ex- 
change diagrams, however, the outgoing A is generated by the handbag part of the diagram. 
If both mesons were on-shell a spin flip of the quarks at each vertex would result in parallel 
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Depolarisation Tensor D N 




FIG. 14. The depolarization tensor Dnn as a function of the Feynman variable xf which measures 
the ration of the actual momentum of the A projected onto beam direction divided by the greatest 
possible one. The error bars on our numerical results represent the statistical error of the Monte Carlo 
integration. 



spins of the incoming proton and the outgoing A. This would be a positive contribution to 
Dnn. Due to the off-shellness of the exchanged meson some small negative contributions 
arise. The actual magnitude depends on the particular kinematical situation considered. 
In essence, the depolarization tensor is controlled by the size of pion and kaon exchange 
contributions and in particular by the phases of the diagrams. These phases are completely 
controlled by the kinematics of the process. We note that in other model calculations these 
phases are either adjusted |38| or interference terms were completely omitted p9| . 

For large negative values of the Feynman-parameter xf the outgoing A is dominantly 
produced by the unpolarized target proton. This causes the depolarization tensor to vanish. 



According to figure TA our model calculations reproduce that feature. For xp > we obtain 
a sizable and positive depolarization tensor. This results from the fact that the leading 
contribution to the process stems form the diagram #1 in figure |8|. Here the pion couples 
to the quark while the diquark acts as spectator. As discussed above such diagrams mainly 
produce A-spins that are parallel to the spin of the incoming proton. On the contrary the 



experimental results suggest that the main contribution should stem from kaon exchange [37 



We note, however, that this obvious discrepancy between theory and experiment has been 
found in other model calculations as well, cf. refs. (B9|, fl40j. 



D. Kaon Photoproduction 

Here we will discuss our numerical results for the process 7^ — > KA. The technical 
details that enter this calculation are given in appendix 0. 

In figure [15] we display the total cross section cr^p — > KA) as a function of the photon 



energy E. We observe that the parameter sets (II-IV) predict cross sections that are com- 
parable with the experimental data. These model calculations do not include axialvector 
diquarks. Once these degrees of freedom are taken into account (sets V and VI), the cross- 
section is overestimated by about a factor four. For the five sets II- VI we find that the 
total cross section is strongly dominated by the kaon-exchange diagram (cf. figure [7]) while 
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the handbag-type diagrams can almost be neglected.0 As only a single diagram contributes 
interference does not occur and it is obvious that the model calculations do not reproduce 
the dip in the energy region l.lGeV < E < 1.4GeV. Tuning the model propagators such 
that the two diagrams are of equal importance this dip could be reproduced p4j. As it has 




1.3 1.5 

E [GeV] 




FIG. 15. The total cross section for kaon photoproduction as a function of the incident photon 
momentum E (left panel) and the differential cross section for kaon photoproduction averaged in the 
energy bin [1.2, 1.3]GeV. The parameter sets (l-VI) are defined in table QTL The experimental data are 
taken from ref. E2|. 



been the case for the associated strangeness production we find that utilizing an exponential 
dressing function (set I) widely overestimates the experimental data. In this case actually 
the handbag diagrams dominate while the kaon exchange contributions are comparatively 
tiny. Figure [15] clearly shows that the large disagreement of the model results with the data 
certainly is not a fine-tuning problem. Rather we must conclude that the comparison with 
data rules out propagators that strongly rise in the timelike region as the one dressed by 
an exponential function does. Contrary to the case considered above the kaon exchange 
diagram exceeds the handbag diagram by almost one order of magnitude for the parameter 
sets II-IV. Possible corrections from subleading meson amplitudes in the kaon-quark vertex 
might decrease the strength of the kaon-nucleon-A form factor in a similar fashion as they 
do for the pion-nucleon form factor. 



In figure [15] we also present the differential cross section in the energy interval 1.2GeV < 
E < 1.3GeV as a function of the angle between the momenta of the initial proton and the 
final kaon in the center of mass frame, cf. appendix [CJ Although the model calculations 
reproduce the empirical increase of the differential cross section as cos 8 cm goes from minus 
to plus unity, the increase appears to be overestimated. For those parameter sets (II-VI) 
for which the resulting cross sections are dominated by the kaon exchange diagrams the 
predicted differential cross sections turn out too small in the backward scattering region 
while they are too big in the opposite direction. As a result the total cross section agrees 
with the empirical data reasonably well. Again, the exponential type propagators yield 



8 It is interesting to note that even for the handbag diagram alone results obtained with a Ball-Chiu 
or bare photon quark vertex, respectively, differ by at most a few percent. 
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differential cross sections that are way off the data and we repeatedly conclude that this 
type of propagators is ruled out. 

We have also computed the asymmetries that are defined in eqs ( |C14| )-( pl6| ). We find 
that they essentially vanish for the model propagators that we consider reasonable, i.e. 
sets 1 1- VI. Although the model calculation correctly predicts that the polarized photon 
asymmetry E, see eq ( |C16|) , is positive for cos 6^ < and negative otherwise, the absolute 
values are off by several orders of magnitude. Only when substituting propagators that are 
characterized by the exponential dressing function the predicted asymmetries roughly agree 
with the empirical data. However, we have discarded already that propagator for other 
reasons given above. 

VI. CONCLUSIONS 

In this paper we have considered baryons as fully relativistic bound states of quarks and 
separable quark-quark correlations, i.e., diquarks. The main purpose of this study has been 
to utilize empirical information in order to restrict the structure of the propagators that 
model confined quarks and diquarks. These propagators enter the four-dimensional Bethe- 
Salpeter equations from which we have computed the mass eigenvalues and wave functions 
that are associated with physical baryons. 

The full covariance of the model wave functions allows us to unambiguously calculate 
form factors up to momentum transfers of several GeV. For spacelike momenta the empirical 
form factors can be very well reproduced with tree-level quark and diquark propagators || . 
On the other hand the description of processes involving timelike momenta is obscured by 
the presence of quark thresholds in the tree-level propagators. Of course, these thresholds 
are unphysical and reflect the absence of confinement. It is thus appropriate to modify 
these tree-level propagators in order to implement the confinement phenomenon. In this 
framework we have considered three qualitatively different cases: In the first case, the tree- 
level poles at timelike real p 2 have been traded for a pair of complex-conjugate poles. In 
this case the imaginary parts (and therefore thresholds) cancel. In the second scenario, the 
pole on the timelike real p 2 axis has been screened at the expense of an essential singularity 
for infinite timelike momenta. In the third case, we have emphasized the issue that the 
propagators should equal the tree-level ones for all complex values of p 2 as \p 2 \ — > oo. 
Together with the condition that no poles occur this property enforces a non-analytic form. 
We have then investigated the phenomenological implications of either of these forms rather 
than attempting to precisely reproduce the experimental data. Obviously, those processes 
are most interesting whose computation involves timelike momenta of the order of one 
GeV entering the model propagators. In diquark models the flavor algebra alleviates the 
calculation of processes with a A hyperon in the final channel. We have therefore focused 
on kaon photoproduction and associated strangeness production with the photoproduction 
being, at least in principle, more sensitive to timelike momentum transfers. 

The model parameters have been fixed by fitting the baryon spectrum and the nucleon 
electromagnetic form factors. Both, the results for the magnetic moments of the nucleon 
and those for the ratio Ge/Gm show that it is important to include contributions from the 
axialvector diquark. The strong form factors g^NNiQ 2 ) and gKNkiQ 2 ) for spacelike momenta 
Q 2 > depend on the amount of admixture of axialvector diquarks in the baryon wave 
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functions. Our numerical result for g n NN(Q) using both, scalar and axialvector diquarks, 
overestimates the empirically determined value by approximately 30%. A possible reason for 
this discrepancy is the omission of subleading amplitudes in the pion-quark vertex. Future 
calculations should therefore include these contributions and the ones from the kaon-quark 
vertex as well. In any event, all these observables are almost insensitive to the specific 
structure of the propagators. Therefore they do not provide an adequate tool to distinguish 
between different parameterizations of the confinement phenomenon. 

The production processes, on the other hand, strongly depend on the form of the prop- 
agators in the timelike region. In particular we have observed that the class of propagators 
that is characterized by an exponential growth for large timelike momenta overestimates the 
cross sections by orders of magnitude. We have associated this failure to the dominance of 
the handbag-type diagram. Apparently any quark propagator that for timelike momenta 
is significantly more enhanced than the tree-level one immediately implies the dominance 
of this diagram. The obvious conclusion is that those propagators should be discarded. 
The other two forms of the propagators have the potential to describe the data reasonably 
well. As mentioned, we have omitted the so far undetermined subleading contributions in 
the kaon exchange diagram that dominates the kaon photoproduction amplitude. From the 
discussion of the pion form factor we have sufficient reason to believe that the inclusion of 
such contributions will favor parameter sets that contain axialvector diquarks. Although 
the non-analytic form for the propagator could not be ruled out by quantitative arguments 
we nevertheless think it should be discarded, because it poses several fundamental problems 
related with gauge and translational invariance which we have detailed in the text. There- 
fore propagators that are characterized by pairs of complex conjugated poles seem to be 
best suited for further studies. 
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APPENDIX A: FROM THE RELATIVISTIC THREE QUARK PROBLEM TO 

THE DIQUARK-QUARK MODEL 

The six-point function G(xi,yi) = (0|Tf|^ =1 q(xi)q(yi)\0) represents the starting point 
for our study of the relativistic three-quark problem. Here the variables Xi and yi not only 
represent the space-time coordinates but also include the discrete labels color, spin, and 
flavor. The six-point function obeys the Dyson equation 



G = G + G ®K®G. (Al) 



The entries of the Dyson equation ( |AT| ) are the disconnected six-point function Go that 
describes the free propagation of three quarks and the three-quark scattering kernel K 
that contains all two- and three-particle irreducible diagrams. The symbol "<S>" in eq (AT) 
denotes summation/integration over all independent internal coordinates and labels. Unless 
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explicitly stated otherwise we will henceforth work in momentum space with Euclidean 
metric. It is thus not necessary to introduce different symbols for momentum and coordinate 
space objects. 

A three-particle bound state with mass M manifests itself as a pole in the six-point 
function at — P 2 = M 2 where P = pi + P2 + P3 is the total four-momentum of the three- 
quark system. We may thus parameterize the six-point function in the vicinity of the pole 

as 

G(khPi) „ t^hMuPiM , (A2) 

where ip denotes the bound state wave-function. Substituting this parameterization into 



the Dyson equation ( |Al| ) and identifying residua, we find the homogeneous bound state 
equation 

ijj = G ®K®ij G- 1 ®ip = 0. (A3) 

Despite its simple appearance this equation is infeasible as neither all two- and three-particle 
graphs, K, nor the fully dressed quark propagator, that is contained in Go, are known. We 
will have to resort to approximations that render the problem tractable and a posteriori 
validate these approximations from the resulting bound state properties. 

The problem greatly simplifies when discarding all three-particle irreducible graphs from 
the interaction kernel K. The kernel may then be written as the sum of three two-quark 
interaction kernels, 

K = K x + K 2 + K 3 . (A4) 

We adopt the notation that the subscript of Ki refers to the spectator quark t&. The respec- 
tive interacting quark pair is (qj, qu) with the three labels (i, j, k) being a cyclic permutation 
of (1,2,3). These two-quark interaction kernels govern the Dyson equation for the two- 
quark correlation functions, go 

g t = G Q + G ®K i ®g i . (A5) 

As the appearance of the free six-point function Go suggests we have defined both and K± 
in the three-quark space. This is easily accomplished by attaching the propagator^ Si of the 
spectator quark to and its inverse S~ l to K^. Expressing eq (|A5|) as Gq®Ki = 1 — Go^g^ 1 
allows us to replace any of the three operators Gq <8> Ki in the bound state equation 



V>= (l-Go®^ 1 ) ® ip + G ® {Kj + K k ) ® tp ip = g i ®{K j + K k )®i) . (A6) 

To further elaborate this form of the bound state equation we define the matrix via 



9 In the framework of these integral equations we factorize the momentum conservation {pi\Si\p'j) = 
(2ir) 4 S(pi)5 4 (p' i — pi). Here S(p) is the ordinary Dirac propagator while Si refers to an operator in 
functional space. We adopt analogous conventions for the other operators. 
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FIG. 16. The Faddeev bound state equation for the component ip\. The equations for ip2 and 
follow by cyclic permutation of the particle indices. 



gi = G + G <g> U <g> G . (A7) 

This reflects the amputation of the external quark legs from after having separated the 
non-interacting contribution. As already mentioned we carry along factors of the quark 
propagators and its inverse to formulate the problem in the three-quark space. For later 
convenience we therefore define 

U = U o S t (A8) 

with the additional factor removed. We have introduced the symbol "o" to denote simple 
multiplications without any contractions because the so-combined operators act on different 
quarks. Finally we introduce the Faddeev components ipi by 

ipi = G ® Ki ® tp . (A9) 



Upon rewriting the definition ( |A.7|) for ti as g^ (g> G = 1 + G <E> U we find the bound state 



equation ip = ipj + ip k + G <8> U <g) (ipj + ipk) an d thus 

ipi = G ® £i ® + V*) = (Sj o 5 fc ) ® ^ ® (V'j + • (A10) 
These are the famous Faddeev bound state equations relating the Faddeev component ipi 



to and ^fc. The graphical representation of these equations is shown in figure [16]. These 
equations embody the full two-quark correlation function ti instead of the kernel K. The 
relativistic Faddeev equations are a set of coupled four-dimensional integral equations and 
represent a considerable simplification over the original eight-dimensional integral equation 
problem defined in eq ( |A3[ ). Unfortunately the Faddeev components ipi still depend on the 
two relative momenta between the three quarks. Expanding these components in Dirac 
space yields an intractable number of coupled integral equations. We therefore further 
simplify the bound state problem. Denoting the incoming and outgoing momenta by re- 
spectively ki and pj we assume that the two-quark correlations ti do not depend on any of 
the scalar products ki ■ Pj that connect momenta of the incoming and outgoing channels. 
This assumption allows us to expand ti in terms of separable correlations in the two-quark 
subspace that is characterized by j,k ^ i and j ^ k, 

t^kj.k^p^pk) = ^2xi(kj,h) D a>a ,(kj + k k ) xt'(Pj,Pk) ■ (All) 

a, a' 

We call these separable correlations "diquarks" and comprise the various types together with 
their discrete quantum numbers within the label a. Note that the propagator D a ^ a i(kj + kk) 
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is diagonal in the discrete quantum numbers a and a' except for the Lorentz indices for the 
axialvector diquark. Furthermore \ a represents the vertex function of two quarks with a 
diquark. Correspondingly x a denotes the conjugate vertex function. The expansion ( All ) 
is pictured in figure [17]. 

In a full solution to the Faddeev problem the will have to be determined from the 
Dyson equation for ti 

ti = Ki + K { <g> G ® U , (A12) 

that follows from eq ( |A~5| ) and involves the kernel components K{. Rather than determin- 
ing these vertices and the diquark propagators from that Dyson equation we will adopt 
phenomenologically motivated parameterizations for these quantities. 

Upon the separability assumption for the two-quark correlations we continue to formu- 
late a relativistic description of baryons based on the Faddeev equations ( |A1(J| ). In this 
approach it is advantageous to introduce an effective vertex function, for the interaction 
of the baryon with the quark and the diquark. This vertex function depends only on the 
momentum of the spectator quark, pi, and the momentum, pj +pk of the diquark quasiparti- 
cle. Eventually this can be reexpressed as a dependence on the relative momentum between 
quark and diquark, pi as well as the total momentum P: 0" = 4>%(pi, P). These dependencies 
are further detailed in refs. H42,43"|. In contrast to the non-relativistic formulation we have 



some freedom in the definition of the relative momentum. We may write 

Pi = Pi - r)P = (1 - v)Pi - viPj + Pk) , (A13) 

where the parameter 77 parameterizes the partition of this momentum among the con- 
stituents. Of course, physical observables like the mass of the bound state or form factors do 
not depend on this parameter (up to numerical uncertainties). The superscript in 0° selects 
a diquark content, a, that builds a baryon together with the quark of species A suitable 
ansatz for the Faddeev components ipi then reads 

1>? 7 = SfsfST' J&v D aa/ <f>( a , , (A14) 

aa' 

where we have made the quark labels explicit. As usual, we sum over doubly appearing 
indices. The quark label i fixes the diquark indices (jk). The momentum routing follows 
these indices as well as the diquark labels a and a' . For further guidance we have visualized 



this ansatz in figure [18]. Noting that Go = Si o Sj o Sk we find the coupled integral equations 



for the effective vertex functions 

E [xl^r'xiya] [DwSfflj,] + (j — k) , (A15) 



i,a 



k ^ ft k J<KJ W> ^ft 

a, a' 



FIG. 17. The separable matrix t\. Also indicated is the amputation of the external legs in the quark 
propagators. 
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FIG. 18. The ansatz ( |A14|) for the Faddeev component ip% of the bound state wave function if) 
using effective baryon-quark-diquark vertex functions <j)\. 



when inserting the ansatz ( |A14j ) together with the diquark parameterization ( |A1 1| ) into the 
Faddeev equations ( |A10|) . In deriving eq ( |A15| ) we have utilized that the quark-diquark 
vertex functions are antisymmetric under the exchange of the quark labels, xtp 1 = ~Xi 1 B- 
This feature is a consequence of the Pauli exclusion principle. We have arranged the terms 
in eq ( |A15| ) such as to exhibit the similarity with the structure of Bethe-Salpeter equations. 
The first term in square brackets represents a six-point function for quarks that is governed 
by the exchange of a single quark between a quark and a diquark. By coupling to the vertex 
function via the propagators for quarks and diquarks it serves as the interaction kernel that 
generates the Bethe-Salpeter equation for a bound state of quarks and diquarks. Thus the 
Bethe-Salpeter equation sums the ladder-type quark exchange diagrams between quarks 
and diquarks. Using the above definitions for total and relative momentum one arrives now 
at the Bethe-Salpeter equation (|6|). 



APPENDIX B: DECOMPOSITION OF THE DIQUARK QUARK 
BETHE SALPETER AMPLITUDE 

Here we will make explicit the full structure of the vertex functions (f>1 for the case of 

the nucleon-quark-diquark vertex. For identical quarks the nucleon-quark-diquark vertex 

functions 0" do not depend on the quark label i. The vertex functions consist of a spinor in 

the case of a scalar diquark (a = 5) and a vector-spinor in the case of an axialvector diquark 

(a = fi — 1 . . . 4). Using positive energy spinors u(P) with P being the nucleon momentum, 

/ $5 \ 

we define matrix-valued vertex functions $ = I ^ J via 

r(p,p) = Hp,p)u(p). (bi) 

Upon attaching quark and diquark legs to $ we obtain the matrix-valued Bethe-Salpeter 

Vf/ 5 



wave functions \T/ 



D( Pd ) := I D( * d) ° ) , (B2) 



D^(p d ) 

*(p, P) = [S(p q ) o D(p d )\ $(p, P) . (B3) 

We demand that the nucleon Faddeev amplitude eq ( |A14j ) has positive parity and de- 
scribes positive-energy states. The latter condition enforces that wave- and vertex functions 
are eigenfunctions of the positive-energy projector A + = |(1 — if/M), i.e., 
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Si A + 

c>2 —if>T-k + 

A^ P^A+ 

A% -iP^f> T K + 



^ Vf(#-Wfr)A+ 



TABLE IV. Basic Dirac components of the nucleon vertex function. The hat denotes normalized 
four-vectors, p = p/\p\. For the complex on-shell nucleon momentum, we define P = P/iM. The 
subscript "T" denotes the transversal component of a vector with respect to the nucleon momentum 
P, e.g. pt = p — (p • P)P. 



$ = $ A + and * = * A + . (B4) 

Using these constraints, the most general structure of $ contains two amplitudes (scalar 
functions) Si(p,P) and S2G0, P) coupling to the scalar correlations and six amplitudes 
A±(p, P), . . . , Aq(p, P) for the axialvector correlations within the nucleon. Explicitly, 

2 

^(p,P) = Y,S l (p 2 ,p-P)S l ( P ,P) , 

i=i 

6 

$"(p, P) = MP 2 ,P ■ P) 7s4- P) ■ (B5) 
i=i 

The Dirac components Si, ■ • • ,A% that obey the positive energy condition are listed in ta- 



ble [TV]. Also, these components have positive parity. We remark that the wave-function \1/ 
can be analogously expanded because it must obey the same constraints as $ does. In the 
nucleon rest frame the individual components of \I/ are eigenf unctions of the three-quark 
spin and orbital angular momentum operators, respectively, when the Faddeev amplitude 
is expanded within the basis ( |B"5| ) [0. Thus, the amplitude Si describes the strength of an 
s-wave in the scalar channel while Ai and A3 represent s-waves in the axialvector channel. 
There also is a small <i-wave component in the nucleon parameterized by A5. All amplitudes 
with even labels are relativistic (lower) components associated with the above described am- 
plitudes that have an odd label, and these flavor components are absent in a nonrelativistic 
description. 



APPENDIX C: TECHNICAL DETAILS FOR CALCULATING PRODUCTION 

PROCESSES 
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FIG. 19. Kinematics for kaon photoproduction pj —* Aif in the center of momentum frame. The 
incoming proton and photon carry the momenta P and p 7 respectively. The outgoing kaon and A are 
labeled by the momenta pk and Pa. 



1. Kaon Photoproduction pj — > AK 

In this appendix we detail the calculation of the diagrams shown in figure [7|. They 
provide the main contributions to the photoproduction process jry — > AK. 

We have performed the calculations in both the rest frame of the proton and in the 
center of momentum system (CMS). For the following discussion we choose the CMS for 
definiteness. The momenta are defined according to figure that is 

P = (-E, 0, 0, iE P ) , p 7 = (E, 0, 0, iE) 
Pk = (\pk\ cos#, \p K \ sin0,O, iE K ) , P\ = (~\pk\ cos6>, -\pk\ sm9,0,iE A ) (CI) 

with 



E P =sjMl + E\ E K = - — {(E + E P f - M\ + M|) 



\p K \ = \J{E K ) 2 -M K , E A =^p K * + Ml. (C2) 

The on-shell conditions and momentum conservation leave only two kinematical variables 
undetermined. These are usually chosen to be the energy E of the incoming photon and the 
angle 9 between the spatial momenta of the photon and the outgoing kaon. 

The three diagrams in figure ^ show the contributions to the transition amplitudes that 
we will discuss here. The uncrossed 'handbag diagram' translates into 

M = if {*a(P/, PA)S( P+ )T K ( q ,p + )S( q )} {r y {p.,q)S(p.)9(p,P)D(p d )} (C3) 

with the momentum routing described in eq (|38D, see also figure |7|. The crossed 'handbag 
diagram' corresponds to the expression 

A 2 = if ^{^A(Pf,PA,)S(p + )T,(q,p + )S(q)}{T K (p^q)S(p^(p,P)D(p d )} . (C4) 
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The definitions for q and p + have changed as compared to the momentum routing for the 
amplitude A\ given in eq (|38l). In eq ( |C4] ) we have instead: 



q = P- - Pk and p + = q + p 7 , (C5) 

with all other momentum definitions unchanged. The amplitude corresponding to the tree 
level diagram arising from kaon exchange is given by 

^3 = (u A [g KNA -f 5 ]u p ) 1 (r%) . (C6) 

q -\- ivi K 

Here Qknk represents the strong form factor that has been discussed in subsection [11 C 




and Mk is the kaon mass. The photon polarization is denoted by e v while T v refers to 
kaon-photon vertex containing the electromagnetic kaon form factor 



T»{q,p) = {qV+p»)F K+ {{p-qy). (C7) 



For the purpose of the present investigation it is sufficient to parameterize the kaon charge 
form factor with a monopole such that the phenomenological value of the kaon radius, 
(r 2 R+ ) = (0.34 ± 0.05)fm 2 , is reproduced: F K+ (Q 2 ) = 1/(1 + Q 2 /(0.84GeV) 2 ). 

In section |Hj we have already reflected on the Bethe-Salpeter formalism in Euclidean 
space. We choose the loop momentum p to be real, thus the temporal components of the 
external momenta in eqs (pl|) - (|Cl]) are purely imaginary. Hence the relative momentum 
Pf as given in (|38"D is complex. Since the solution of the Bethe-Salpeter equation provides 
the vertex function $a only for real relative momenta Pf, we have to extrapolate $a to 
complex momenta. We fit rational functions to the vertex functions that are known at N 
real meshpoints. These rational functions can then easily be analytically continued. For 
real momenta a comparison of the fitted parameterization to the known results allows us to 
estimate the reliability of this treatment. 

The differential cross section depends only on the energy E of the incoming photon 
and the angle 9 between the spatial photon and kaon momenta. That is illustrated in 
figure [19|. The differential cross section is defined with respect to the solid angle element 
dflx = 2ird(cosd) of the outgoing kaon: 

da , „ „ s 1 ^ — - da 



dflx 4 ^— ' dVtx 

with 

da ,„ „ s a 



, o (£,#) = __ AM P M A \A 1 + A 2 + A 3 \ 2 . (C9) 

We average, respectively sum over the spins Sj, Sf of the initial and final states. The phase 
space factors denoted by a are given as 



1 \p K \ 2 d\p 



P- Py E K E A d(E K + E A ) (C10) 
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with 



d\p 



K 



d(E 



K 



\p 

E 



K 



K 



Pk ■ Pa 
E^ 



-i 



(Gil) 



Note that the right hand side of eq (|C11| ) remains positive given that > Mr- In the 
CMS the expression ( |U1UD for a simplifies to 



a 



Wk\ 1 



(C12) 



where s = (P + p 7 ) 2 = (j>k + -Pa) 2 denotes the total center of mass energy squared. In 
obtaining the phase space factors (|C10| ) we adopted the one particle normalization conditions 
(for Minkowski space) 



p u 



(p\p') B = 2p°{2n) 3 5 3 {p- p>) and (p\p') F = — (27r) 3 5 a (p - p') 



(C13) 



for bosons (B) and fermions (F). These conventions also enter the calculation of the transi- 
tion amplitudes z = 1,2,3 and the normalization of the Bethe-Salpeter wave-functions. 

We obtain the various asymmetries by restricting the sum in eq flCB|) over the spins to 
two of the three non-scalar particles. We thus obtain the A-polarization asymmetry 



P(E,9) = - 
the polarized target asymmetry 

and the polarized photon asymmetry 



E(E, 



1 



E 



E 



E 





=T] - 


[sa 


=1] 


[sa 


=T] + 


[sa 


=1] 




=tl- 


[s p 


=1] 


1 S P 


=T] + 


[Sp 


=1] 


[s 7 - 


=tl- 


s 7 = 


=1] 


[s 7 = 


=t] + 


> 7 = 


=1] 



(C14) 



(C15) 



(C16) 



where we used the shorthand notation 



[sa =T] 



da 



dtt 



K 



etc. 



(C17) 



SA = T 



Furthermore we denote the spins of the photon, the proton and the A by s 7 , s p and sa, 
respectively. 

The total cross section is finally obtained from eq (ICS) via 



a{E) 



dtt 



K 



(C18) 
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FIG. 20. Kinematical variables for the reaction pp — > pKA in the center of momentum system. 



2. Associated Strangeness Production pp — * pKA 

Again we calculate the cross section for the productions process pp — > pA'A in the center 
of momentum frame. The kinematical setup is depicted in figure (|20D and amounts to the 
momentum routing 

P x = (|pl cos(^ a ), \p\ sin(0 A ), 0, iE) , P 2 = (-\p\ cos(# A ), -\p\ sin(0 A ), 0, iE) , 
if = (i^, \K\ sm(9 K ) cos(4> K ), \K\ sm(9 K ) sin((f> K ), iE K ) , 

A = (A x , 0, 0, iE\) , Pf={P f ,iE P ), (C19) 
where i£ denotes the center of momentum energy. We have furthermore introduced 



E P = 2E-E*-E K , A X = JE 2 A -Ml, \K\ = J E\ - M\ , P f = —K — A 



El-El-kl + Ml-Ml 
&x = ^ , sm{V K ) 



\ 




(C20) 



where all the momenta are Euclidean. There are now four independent variables left: E\, 
Ex-, 0\ and 4>k- As for pr) — > KA we find the vertex function r A for complex arguments by 
analyticly continuing a fitted rational function. 
The differential cross section is given by 

^-^E l <m\M s , s ,, r , r r\ 2 , where rfS = — *- , M } M ^ dE K dE K d<\> K , (C21) 
dS7 A 4 128vr 5 p cm \E 

with fif2 A = 27r dcos(^A)- The masses of the proton and the A are denoted by Mp and M\ 
while (p cm ,iE) represents the four vector of the total momentum in the CMS. The ampli- 
tude Ai s ,s',r,r' f° r the reaction py — > KA depends on the spin orientation of the incoming 
and outgoing particles. In ( |C21| ) we average, respectively sum over the spin projections of 
the incoming and outgoing particles. The integrations over the kaon and A energies are 
constrained by the available energy, which is 2E. 
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For each diagram in figure |8| and for each specific diquark content (cf. Table |) one has 
a contribution of the form 

c (orb?) H (C22) 

to the amplitude Ai. Here refers to the mass of the intermediate pseudoscalar meson, 
i.e. m<f> = M n or m$ = Mr. The 'meson matrix element' is given by 

C s , s > = u s '(P 2 ) i^g-KNNiQ 2 ) u s (P f ) or C s ,s' = u s >(P 2 ) H5 9kna(Q 2 ) u s (Pa) , (C23) 

depending on whether the intermediate meson is a pion or a kaon. The 'handbag part' 7i 
has the general structure 

n = i J -0- {*A(p/,Pf) S( P+ )T K (p + ,q) S(q)} {r w (g,p_) S(p„)D(p d ) <f> P ( Pi , P^} . (C24) 

Here the Bethe-Salpeter amplitudes and $p for the A and the proton as well as the 
(di) quark propagators S and D enter. The meson-quark vertices 1^ and Tk are defined in 
eq©. 

The calculation of the amplitudes involves four dimensional loop integrations. All the 
momenta are Euclidean and therefore we use hyperspherical coordinates. For the inner 
loop of the handbag diagrams we use a Gauss-Legendre routine, whereas the phase space 
integration is performed with Monte-Carlo methods. Due to the considerable effort it takes 
to integrate eight integrals numerically the calculation could only be performed to an overall 
accuracy of 5% to 15%. However, we consider that sufficient for a comparison with data. 

The depolarization tensor is defined as 



D NN (x F ) = —— (C25) 
a + o 



with the shorthand notation 



« = ^(Ua + Ua) and b=^-(-\ p U + i p U) (C26) 

dXF UXf 

for the cross sections with different spin projections. Here xf denotes the real momentum 
of the A scaled by the maximum value allowed by the kinematics: 

Xf = ] A|11 . (C27) 

-A- 1 1 m.n.T, 
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